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Abstract 



A coarse-grained cellular automaton is pro- 
posed to simulate traffic systems. There, cells 
represent road sections. A cell can be in two 
states: jammed or passable. Numerical cal- 
culations are performed for a piece of square 
lattice with open boundary conditions, for 
the same piece with some cells removed and 
for a map of a small city. The results indi- 
cate the presence of a phase transition in the 
parameter space, between two macroscopic 
phases: passable and jammed. The results 
are supplemented by exact calculations of the 
stationary probabilities of states for the re- 
lated Kripke structure constructed for the 
traffic system. There, the symmetry-based 
reduction of the state space allows to par- 
tially reduce the computational limitations of 
the numerical method. 



1 Introduction 

In simulations of traffic systems, cellular au- 
tomata (CA) are a common tool [H El [3l H] . 
A cellular automaton consists a structure 
of cells, a set of cell states and the rule of 
time evolution which transfers a state of a 
cell with its neighborhood to the state of 
this cell in the subsequent time moment. 
In this description, states, space and time 
are discrete, while the traffic systems at 
lest space and time are inherently contin- 
uous. Still, there is a rich variety of CA 
which enable to investigate properties and 
phenomena of traffic systems, reproduced 
sometimes with surprisingly subtle details. 
As a rule, the traffic CA are classified as 
single-cell or multi-cell models, where a 
vehicle occupies a single cell or more cells jl] . 
On the other hand, traffic networks are also 
parametrized in different ways, as a node 
can represent a stop, a cross-road or a route 
[5]. So simplified when compared with real 



systems, the technology of CA suffers known 
computational limitations: a more detailed 
description is paid by the smaller size of a 
simulated system. 

More recently, a modification of CA has 
been applied where a cell represents a state 
of the whole considered system [6l [7]. This 
approach can be seen as an example of the 
concept of Kripke structures [H] , where nodes 
represent states of the whole system and 
links represent processes leading from one 
state to another. The obvious drawback 
of this parametrization is that the number 
of nodes increases exponentially with the 
system size. In p], this difficulty is evaded 
by taking into account only the states which 
appear during the time evolution. The 
same idea was developed into a technique 
of time series analysis, known as recurrence 
networks. Briefly, the rule of time evolution 
is used to generate new states which are 
attached as nodes to the simulated network; 
in this way the signal is characterized in 
terms of a growing network. For details see 
[0] and references cited therein. 

Our aim here is to discuss a new cellular 
automaton designed for modeling jams in 
traffic systems. The novelty of this automa- 
ton is that cells represent sections of road 
which can be either jammed or passable. 
A jam can grow at its end and flush at its 
front; the competition between these two 
processes depends on the local topology of 
the traffic network. Our description, inspired 
by percolation, is more coarse-grained, than 
in other models. According to the classifi- 



cation of traffic models, presented in [TU] . 
our model belongs to macroscopic queueing 
models. Some model elements remind the 
cell transmission model by Daganzo [HI |12] : 
namely, the rates of inflow and outflow in 
the cell transmission model are similar to 
the rates of grow and flush of traffic jam, 
defined below. However, as it is explained in 
details in the next section, it is only jammed 
and passable cells what is differentiated here, 
and the flows of vehicles are not identifled. 
The price paid is that a range of dynamic 
phenomena as synchronization and density 
waves are excluded from the modeling. 
These phenomena, essential at the scale of a 
road [131 HH Uni [IS], can be less important 
at the scale of a city. Consequently, our 
approach should be suitable for macroscopic 
modeling of large traffic systems. 

Our second aim is to construct the Kripke 
structure on the basis of the same cellular 
automaton. This in turn limits again the size 
of the system, because of the exponential 
size dependence of the number of states. We 
are going to demonstrate that our recent 
tool, i.e. the symmetry-induced reduction 
of the network of states [TTl HHj, is useful 
to partially reduce the computational barrier. 

In the next section we describe the automa- 
ton in general terms and we recapitulate the 
method of reduction of the network of states, 
mentioned above. As the direct form of the 
automaton rules depends on the traffic sys- 
tem under considerations, the exact descrip- 
tion of the rules is given in Section 3, together 
with the information on the analyzed sys- 
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terns, both artificial (the square lattice) and 
real (a small city). Two last sections are de- 
voted to the numerical results and their dis- 
cussion. 

2 The model 

2.1 Cellular automaton 

We analyze a simple automaton, where each 
cell i- a road section - can be either in the 
state s = or s = 1. The state = means 
that a fluent motion via a given road section 
is possible, while the state s = 1 means a 
traffic jam. As each road section is a part of 
larger system, the cell state depends on the 
state of roads where one can enter from a 
given road section. Namely, the probability 
of a traffic jam back propagation as well as 
the probability of a traffic jam to be flushed 
depend on the number and state of neighbor- 
ing road sections. To initialize calculations, 
one has to assign values of three parameters. 
Two of them, w and v, describe the whole 
system, and the last one p is related to the 
boundaries. Specifically, w is the probability 
that a traffic jam arises on a given road 
section due to its presence on the roads 
directly preceding the currently considered 
one (jam behind jam), v is the probability 
of a jam flush (jam behind passable gets 
passable), and p is the probability that a 
traffic jam appears at a road section at the 
boundary, but out of the system. The latter 
parameter describes the system interaction 
with the outer world. The parameters w 
and V can be related to the flows used in 



[T^ I2U] for the discussion for congestion near 
on-ramps. 

The probability of a change of the state 
of a given road section is obtained as the 
result of the analysis of the state of this 
section and the state of its neighborhood. 
We ask for which ranges of the parameters 
the system is passable in the stationary state. 

The detailed realization of the model de- 
pends of the topology of the traffic network. 
In Section 3 the exact algorithm is presented 
together with the presentation of the ana- 
lyzed systems. 

2.2 The state space and its re- 
duction 

The automaton defined above can be used 
for simulations, and the results of these sim- 
ulations are reported below. The same au- 
tomaton is used here also to construct the 
network of states, as in |6l |T7l [18]. This 
network, equivalent to the Kripke structure 
[8|, is formed by all possible combinations of 
states of roads which play the role of nodes. 
Next, an appropriate master equation [21] 
is constructed, which reflects all possibilities 
of states obtained from the current state. 
The obtained matrix of transitions between 
states, i.e. the transfer matrix, is equivalent 
to the connectivity matrix of our state net- 
work. For a given set of model parameters 
{p, w and v), eigenvector of the matrix asso- 
ciated with the eigenvalue equal 1 serves to 
calculate probabilities of particular states in 
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the stationary state. Having these values, one 
can evaluate how passable the system is un- 
der given conditions from the average number 
of unjammed (s = 0) states 

2iv 
i=l 

where: is the size of the system (the 
number of considered road sections), Pi - 
probability of i-th state and nj(0) - number 
of zeros in i-th state. We note that in this 
equation, me make an average over the states 
of the whole network, and not over the states 
of local cells. 

As the obtained number of states is large 
even for moderate systems, we reduce the sys- 
tem size by the application of the procedure 
proposed in our previous papers [T7l[18]. The 
method of the reduction of the system size is 
based on the symmetry observed in the sys- 
tem, which manifests in the fact that prop- 
erties of some elements of the system are ex- 
actly the same. The starting point is the net- 
work of states, and the core of the method 
is to divide nodes into classes; the stationary 
probability of each node in the same class is 
the same [T7l|18]. To begin, for each node the 
list of its neighboring nodes is specified, with 
the consideration of weights of particular con- 
nections. Provisionally, the class of each state 
is determined by its degree; for each state 
its symbol is replaced by the symbol of class, 
which discriminate nodes which have differ- 
ent number of neighbors. At the next stage 
we examine the lists of neighboring nodes in 
terms of class symbols assigned to a particu- 



lar neighbors and weights of appropriate ties. 
If for nodes assigned with the same symbol 
the symbols assigned to its neighbors are dif- 
ferent or their are the same but their weights 
are different, an additional class distinction is 
introduced. At the end of the algorithm, the 
classes, i.e. subsets of nodes are indicated, 
which have identical lists of neighbors with 
respect of the number of neighbors, the sym- 
bol assigned to each of them, and weights of 
particular connections [TTf \TE\. 

3 Analyzed systems 

3.1 Square lattice 

As a reference system we analyze a system 
of directed roads placed on edges of a regu- 
lar square lattice. The lattice is finite, with 
open boundary conditions. For such a sys- 
tem each road has two in-neighbors and two 
out-neighbors. As it will be explained in de- 
tail, the probability of the state change de- 
pends only on the state of out-neighbors. At 
the boundaries, a road has one or none out- 
neighbors (none at the corner). At each road, 
the traffic takes place only in one direction, 
say upwards and right. This setup is bor- 
rowed from the Biham-Middleton-Levine au- 
tomaton [22]. The algorithm of a change of 
the state of the road for the square lattice is 
presented in Figjl} 

In the above algorithm s is a state of a 
given road, Ln is the number of roads given 
road is connected to, the quantity s„ refers 
to the state of the road which is a neighbor 
of the currently considered one (if the road 
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if s = then 
if Ln = then 

else if Ln = 1 then 

if s„ = then 

w 

^ 1 

else 

w w 
P=p— + — 
9 9 
^ — ^ 1 

end if 

else 

if s„i ^ s„2 then 
p= — 

^ 1 

else if s„i = s„2 ~ 1 then 

0- ^1 
end if 

end if 
else 

if Ln = then 

else if L„ = 1 then 
if Sn = then 

V 

P=-{2-p) 

1 ^ > 

else 

V 

1 — ^ ^ 

end if 
else 

if Sni ^ Sn2 then 

1-^0 
else if Sni = Sn2 = 1 then 

V 

p=- 

1 ^ 

end if 
end if 
end if 



has two neighbors their states are marked 
respectively as Sni and s„2)- The probabihty 
P of a change of the state depends, as 
it was mentioned above, on the state of 
the considered road and the state of the 
neighborhood which determine how probable 
the change is. Namely, the transition from 
jammed to passable (1 to 0) mean that the 
jam at a given road section is flushed by free 
motion of vehicles at the jam front. This is 
possible only if the out-neighboring section is 
empty. Further, the transition from passable 
to jammed (0 to 1) is possible only if the 
out-neighboring section is jammed. In both 
cases, the transition depends on the state of 
the out-neighbors; the state of in-neighbors 
is not relevant. 

In Fig|2] a piece of the system is presented. 
Each road can be either passable (in our no- 
tation a road is in a state 0) which is marked 
as a dashed-line or a traffic jam can be formed 
(a road is in a state 1) which is marked as a 
solid-line. The direction of traffic was ticked 
on roads in the state 0, but the rule is the 
same for all roads; we keep down-up and left- 
right direction. Here we present one of the 
possible changes of the state of the system. 



Figure 1: Algorithm of a change of the state 
of the road for the square lattice. 



3.2 Small city 

The method was also applied to a real road 
network of a small Polish town Rabka. The 
structure of roads which matter in traffic was 
selected - dead ends are removed (Figj3|. 
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Figure 2: Example of the change of traffic. A 
dashed-hne refers to the state (fluent flow), 
and a solid-line refers to the state 1 (traffic 
jam). Arrows indicate the direction of traffic. 



Each road, if necessary, was divided into sec- 
tions of approximately equal length. We get 
374 sections. Here the number of neighbors 
for different roads varies as it results from 
the town topology. Each section is a two-way 
street. In this case the algorithm has a form 
presented in Figjij 

In the algorithm presented in Figjl] sum- 
ming goes through the states of the roads out- 
going from a given one, and kout is a number 
of outgoing roads. 

4 Results 

4.1 The square lattice 

All presented results are a time average in 
the steady state over 100 realizations for the 
square lattice of the size = 100 x 100. To 
check that the results do not depend on the 
initial conditions, we use three options for the 
initial state: states of all roads set to 0, states 
of all roads set to 1, and a state of each road 
is set randomly to be or 1. 
The results depend on the values of the model 



Figure 3: Road network of Rabka. Roads 
are divided into 374 approximately equal sec- 
tions. The black dots mark the exit /entrance 
roads. 



if s = then 

if Ln = then 

w 



else 



'l^^l^ 1 

end if 

else 

if Ln — then 

else 

P=vJ:(1-s) 

1 > 

end if 

end if 

Figure 4: Algorithm of a change of the state 
of the road for a real network. 
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parameters p, w and v. As a result for the 
whole system the percentage of roads in the 
state (#0[%]) is calculated. The higher the 
number of zeros the more passable the system 
is. The results for two different values of the 
parameter p are presented in Fig. [5] for p = 
0.1 and for p = 0.7. 



interpreted as an indication of a phase tran- 
sition. To verify its sharpness dependence 
on the system size, we calculated the curve 
#0 vs V for a selected case: p = 0.7, w=0.5 
and various system sizes N"^. The results 
are shown in Fig|6} Indeed, the sharpness 
increases with A^, and the curve for = 200 
is close to the step function. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



(a) p = 0.1 




(b) p = 0.7 

Figure 5: Diagram for the square lattice of 
the size 100 x 100 (average in the steady state 
over 100 realizations). 



The increase of the percentage of zeros 




3.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0. 



Figure 6: Percentage of zeros in a function of 
V for square lattices of different sizes N x N 
for p = 0.7 and w = 0.5 (average in the steady 
state over 100 realisations). 



We also check how removal of some num- 
ber of roads changes the obtained results, to 
check if the symmetry of the square lattice 
is necessary. In FigjT] the results obtained 
when 10% randomly chosen road sections is 
removed. The removal is done separately for 
each realization. If, in a consequence of the 
removal, some part of lattice is isolated, it is 
removed as well. The results, shown in FigjT} 
indicate that the phase transition, found for 
the square lattice, is observed also in a ran- 
domized lattice. The maximal number of ze- 



with the parameter v, visible in FigjS} can be ros in this case is less than 90 percent, be- 
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cause the plot is normalized to the whole 
square lattice, including the removed links. 




(a) p = 0.1 



The main difference between this network, 
as constructed from the map in FigjSj and 
the square lattice (with removals or not) is 
that the Rabka network is less connected. 
There, often the road sections form long 
chains. Comparing FigsjS] and [7] we see 
that the consequence of this difference 
is that jammed state is less likely. The 
origin of this result is that jams are created 
behind the jammed road sections; the more 
in-neighbors of these sections, the more 
jams appear. Besides of that, the obtained 
plot are similar to those for the square lattice. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



(b) p = 0.7 

Figure 7: Diagram for the square lattice of 
the size 100 x 100 (average in the steady state 
over 100 realisations), with removal of 10% 
randomly chosen road sections. 



4.2 Small city 

The result obtained for the simulations of 
the traffic network in Rabka, formed by 
374 road sections, are presented in FigsM 



4.3 A simplified map 

Exact calculations of P{0) can be performed 
merely for systems much smaller than 
hundreds of road sections. For the sake of 
comparison of the methods, we simplified 
the map of Rabka leaving only nine two-way 
roads. This leads to the system of 2-'^'^ states. 
The results of the simulations for this system 
are shown in [9j As the system is much 
simplified, the results for the full (Figjsj) 
and reduced (Fig|9]) traffic networks differ 
substantially for p = 0.1. Surprisingly, those 
for p = 0.7 are quite similar. The same 
simplified network is solved exactly by the 
solution of 2^* = 262144 Master equations 
pi] for the stationary state for different 
sets of the model parameters of the related 
Kripke structure. In this exact method, the 
parameters p, v, w enter to the weights of 
links between states, or - equivalently - to the 
rates of the processes which drive the system 
from one state to another. For each case 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 



(a) p = 0.1 



(a) p = 0.1 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



(b) p = 0.7 



(b) p = 0.7 



Figure 8: Diagram for the real network of the Figure 9: Diagram for the real reduced net- 
size N = 374. work of the size N = 18. 
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we can then calculate, in accordance with 
EqjT| the mean stationary probability that 
the road sections are passable. Obtained 
results are presented in FigsfTOj The same 



figure shows the solution for the classes of 
states, as described in [TTl IB]- In this case, 
the class identification procedure allows for 
the reduction of the system size about twice, 
to 102400 classes. 



5 Discussion 

The goal of this paper is to describe large 
traffic networks with a cellular automaton, 
where states of road sections are reduced 
to two: passable and jammed. The coarse- 
grained character of the new automaton is 
close in spirit to the percolation effect. The 
results of our simulations allow to identify 
a phase transition between two macroscopic 
phases, again passable and jammed. Addi- 
tionally, the calculations are repeated for a 
much smaller traffic network, constructed 
by a strong simplification of a map of a 
small Polish city. These calculations are 
performed to compare the results with 
the exact solution of the stationary state, 
obtained by two equivalent methods. This 
comparison suggests, that the accordance of 
simulation with the exact solution is better 
for more jammed systems, i.e. more close to 
the phase transition. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



(a) p ^ 0.1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
V 

(b) p = 0.7 

Figure 10: Diagram for the real reduced net- 
work of the size = 18, obtained from the 
analysis of the state space. Here cl refers to 
the network of classes, and st to the network 
of states. 



The drawback of our automaton is that all 
information about specific local conditions 
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of traffic jams cannot be reproduced. Tfie 
model captures merely the jam spreading. 
The parameters w and p depend on the 
external state, and serve as an input for the 
calculations. The parameter v should be 
cahbrated separately for each traffic system. 
After this calibration, the main result of 
the model - the probability of the jammed 
phase - should be reproducible and useful to 
control the traffic phases. The advantage of 
the model is its simplicity, which allows to 
to simulate larger traffic systems in real time. 
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